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Abstract 

Dendritic growth, and the formation of material microstructure in general, necessarily involves 
a wide range of length scales from the atomic up to sample dimensions. The phase field approach 
of Langer, enhanced by optimal asymptotic methods and adaptive mesh refinement, copes with 
this range of scales, and provides an effective way to move phase boundaries. However, it fails to 
preserve memory of the underlying crystallographic anisotropy, and thus is ill-suited for problems 
involving defects or elasticity. The phase field crystal (PFC) equation — a conserving analogue of 
the Hohenberg-Swift equation — is a phase field equation with periodic solutions that represent 
the atomic density. It can natively model elasticity, the formation of solid phases, and accurately 
reproduces the nonequilibrium dynamics of phase transitions in real materials. However, the PFC 
models matter at the atomic scale, rendering it unsuitable for coping with the range of length scales 
in problems of serious interest. Here, we show that a computationally-efficient multiscale approach 
to the PFC can be developed systematically by using the renormalization group or equivalent 
techniques to derive appropriate coarse-grained coupled phase and amplitude equations, which are 
suitable for solution by adaptive mesh refinement algorithms. 

PACS numbers: 81.16.Rf, 05.10.Cc, 61.72.Cc, 81.15.Aa 



I. INTRODUCTION 



During the last thirty years or so, the field of computational materials science has emerged 
as a flourishing sub-discipline of condensed matter physics. It is now relatively straight- 
forward to compute realistic-looking materials microstructures for a variety of processing 
conditions, and, with enough computer power, to begin to make quantitative predictions 
about phase diagrams, morphological phase diagrams, growth rates and other probes of the 
kinetics of phase transitions. But this was not so, nor even obviously within the realms 
of possibility, during the 1970's when Jim Langer, Pierre Hohenberg and others initiated 
studies of phase transitions kinetics and instabilities in spatially-extended systems, whose 
fruitful union in the 1980 's gave birth to the field of pattern formation as practiced today. 
In those early days, the key problems were to understand the nature of the instabilities 
around uniform or similarity solution states, and to capture correctly the intricate non-local 
and retarded feedback between interface dynamics and the (typically) diffusion fields around 
them. To paraphrase John Archibald Wheeler's epigram on general relativity: solidification 
fronts tell heat how to flow; heat tells solidification fronts how to move. 

Following a variety of innovations from Langer and collaborators, especially the in- 
troduction of phase field modelsQ, Q], and their practical implementation with improved 
asymptotics S] and adaptive mesh refinement 4^ these difficulties are largely resolved, cer- 
tainly in principle, and to a significant extent in practice 5] . Three dimensional structures, 
including fluid flow effects, can now be simulated on desktop computers 0, Q|. The key con- 
ceptual difficulty that remains as a challenge to this day is the huge range of length and 
time scales encountered in pattern formation processes. For example, in solidification, one 
encounters the capillary length on the scale of 10~ 8 m and the dendrite tip radius on a scale 
ranging from tens of microns to millimetres, depending on the undercooling. The latent heat 
around a solidification front extends to a distance known as the diffusion length, which can 
be 10 _4 m or larger; and finally, in a real processing experiment, there will be the system 
size itself, many orders of magnitude larger still. 

The problem of "bridging the length and time scales" between atomic and sample di- 
mensions is the focus of the present article. Despite much activity to address the scale- 
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variants}^, most existing work is currently limited to crystalline materials with a few iso- 

as an exception a promising approach [22] based on the phase 
field approach, due originally to Langer (but not published for many years) [l]. 

The approach advocated here introduces the atomic level of description into phase field 
approaches using a recently developed formalism known as the phase field crystal (PFC) (23! 



241 ]. The PFC describes atom crystalline structure as a periodic density wave, and posits a 



natural equation of motion for the density field. Having the character of a continuum partial 



differential equation, it can be coarse-grained using renormalization group (RG) 
related methods (see, e.g. |27ip . developed for the quite different problem a 
hydrodynamic instabilities 28] in spatially-extended dynamical systems 
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The coarse-grained counterparts of the density field are its amplitude and 
phase, for which we derive equations of motion. Once known, the density field can then 
be reconstructed. In this article, we sketch the derivation of effective equations at the 
mesoscale, and show that the numerical solution of the renormalization group equations 
generates solutions that are virtually indistinguishable from brute force solutions of the 
PFC equations of motion. 

A further aspect of our work is that computational efficiency is in practice best achieved by 
using adaptive mesh refinement. The basic idea is that for any field whose spatial variation 
is essentially uniform everywhere, with localized regions of rapid variation, it makes no 
sense to use a uniform mesh in numerical calculations. A coarse mesh can be used in places 
where the field is spatially uniform, and a coarser one is used in the transition zones. In 
the context of phase field models, we have developed techniques to implement these ideas 
in an efficient manner, and demonstrated that the computational complexity of this sort of 
algorithm scales not with the system's volume, but with the surface area of solidification 
front (transition zones) — a considerable saving for large systems J] . The renormalization 
group equations derived below have solutions with the desired character: thus we intend in 
future work to develop adaptive mesh refinement codes for solving these equations. This 
topic is briefly discussed at the end of this article. 
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II. THE PHASE FIELD CRYSTAL MODEL 



The phase field crystal model for a sin gle component system describes the space-time 



behavior of the density p(x, t) and is capable [23|,|24j of capturing realistic aspects of materials 



dynamics, including grain growth, ductile fracture, epitaxial growth, solidification processes, 
and reconstructive phase transitions. By construction, the stationary states of p(x) are 
periodic, and distortions of the density field by external perturbations applied to boundaries 
or by defects, for example, result in a raising of the free energy in accord with Hooke's law, 
with higher order terms representing non-linear elasticity. This feature makes the PFC an 
ideal tool with which to explore nanoscale strain effects and their influence upscale to the 
continuum. 

Let F{p} denote the coarse-grained free energy functional whose minima correspond to 
the equilibrium (lattice) state of a (^-dimensional system, and whose corresponding chem- 
ical potential gradient drives the dynamics of p. For single component, two-dimensional 
systems with a hexagonal lattice state, the appropriate form of F is that originally due to 
Brazovskiij^]: 

F{p(x)} = J d d x [p (aAT + A (q 2 + V 2 ) 2 ) p/2 + V/4] (1) 

where a, A, q a and u can be related to material properties^. I^ . and AT denotes the 
temperature difference from some reference high temperature. For convenience it is useful 
to rewrite this free energy in dimensionless units, i.e., x = rq a , ip = p\juj Xq%, r = aAT/Xq^ 
and F —>■ Fuj X 2 q^~ d . We also use dimensionless time units t, related to the physical time t 
in terms of the phenomenological mobility r by t = TXq^t. In these units, the equation of 
continuity for the density becomes 

d^/dt = V 2 ( [r + (1 + V 2 ) 2 ] t/j + + C (2) 

The conserved Gaussian noise, required by the fluctuation-dissipation theorem to satisfy 
(((ri, ti)((f 2 , r 2 )) = — £V 2 5{fi — r 2 )5(r 1 — r 2 ) with E = uk B Tq li ~ 4: /X 2 , will not generally be 
important for describing phase transition kinetics, and so will henceforth be neglected here. 

Elder and Grant [2^ have studied the mean field phase diagram of the PFC equation 
@ in the one mode approximation, that is valid in the limit of small r, and represented 
in the plane of dimensionless temperature, r, and dimensionless average density, ip. There 
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are three possible equilibrium solutions: a 'liquid', ipc — a two-dimensional 'crystal' 
with triangular symmetry, ipT = At (cos^qrx) cos(grZ// \/3) — cos(2gj , l// \/(3))/2^ +-0, and a 
smectic phase which will be ignored for present purposes. The triangular lattice can exhibit 
defect structures during the relaxation to equilibrium, thus capturing the kinetics of phase 
transformations. 



III. COARSE-GRAINING OF THE PHASE FIELD CRYSTAL EQUATIONS 

The dynamics of the slowly-varying amplitude and phase describes fluctuations about 
a given set of lattice vectors, but must be covariant with respect to rotations of those 
lattice vectors, so that we can properly describe polycrystalline materials with arbitrarily 
oriented grains. A similar situation arises in describing amplitude and phase variations of 
convection rolls, and in the context of the model Swift-Hohenberg [2^ equations, the form 
of the governing equations was originally proposed by Gunaratne et al. J3|, and derived 
systematically from the RG formalism of Chen et al.j^ by Graham j^J (see also ref. jsj]). 

The details of the corresponding calculations for the PFC will be given elsewhere, but 
the results can be quickly obtained by considering the linearized equation governing small 
variations about the constant solution -0: 



dt 



{r + (l + V 2 ) 2 } + 3^ 2 ] ^ 



(3) 



We consider a one- mode triangular perturbation with wavenumber q t , given by ip = 
exp(ut) cos(q t x) cos(q t y/\/3) — cos(2q t y / \/3) /2 . Substitution into Q yields the following 
dispersion relation: 

1 . r / 1 \ -' 

(4) 



r + 3^ 2 + ( 1 - 



The fastest growing mode has q t = i-y/3/2 for this perturbation, which, as was shown by 
Elder and Grant |2J| working in the one- mode approximation, is the mode which minimizes 
the free energy functional Thus, the uniform phase becomes unstable to hexagonal 

perturbations when \tjj\ < J—r/3. 

The triangular phase solution can be written in more general form as 



ijj(x) = Aj(t) exp(ikj ■ x) + ip, 



(5) 
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where k\ = k (— iy/3/2 — j/2), k 2 = k ] and k 3 = fc (iv^3/2 — j/2) are the reciprocal lattice 
vectors, i and j are unit vectors in the x- and y-directions, k is given by q t and Aj(t) = 
Ajexp^-t). 

Substituting into (JBJ) and applying orthogonality conditions, we readily obtain 

^• = -N 2 [-r + {i-N 2 } 2 ], (6) 

where V = — (r + 3ip 2 ), and we have transformed to units where ko = 1. In order to obtain 
the spatial component of the amplitude equation we study the effect of applying small spatial 
modulations to the amplitude, i.e. 

Aj{t) i — ► Aj(x, y, t) = A°j exp(w i (Q)*) exp(iQ • x), (7) 

where Q = Q x i + Q y j is the perturbation vector. The exponent controlling the growth 
rate along each lattice basis vector is given by 

2" 



r-{i- IQ + k^y 



ojj{Q) = \Q + k,\ 2 
Replacing Fourier variables by gradient operators, we obtain 

\Q + kj\ 2 = -V 2 - 2ikj ■ V + 1. 



(8) 



(9) 



Putting together Eq. (JHj) and Eq. (0), we find that space-time amplitude variations along 
each basis are governed by: 



dt 



1 - V 2 - 2ikj ■ V 



r - {V 2 + 2ikj ■ vf 



A r 



(10) 



The manifestly rotationally covariant operator on the right hand side of (jl(Jj) will be denoted 
by Cj. 

The non-linear component of the amplitude equation is obtained bv renormalizing the 
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3 it ] . Expanding 



secular terms in the perturbation series about the periodic state 
about the one-mode triangular phase solution, scaling ip by the factor y/—r, and renaming 
the new variable ip, yields the PFC equation in the form 

dip 



dt 



V 2 (l + V 2 ) 2 ^ + eV 2 (^-^) 



(11) 



where e = — r. Expanding ip as ip = ipo + + ^"^2 + C( e3 ) an d substituting in Eq. (JHJ), 
using Eq. (JHJ) for if) Q , we obtain at 0(e), 



dipi 
dt 



V 2 (l + V 2 )> 1 = V 2 (^-^o). 



(12) 



The secular terms on the right hand side can be identified from the expansion 



V 2 (V 3 - Vo) = (1 - 3^ 2 ) E A i ex P(^ • 

j 

-3A 1 (\Ai\ 2 + 2\A 2 \ 2 + 2|A 3 | 2 ) exp(ijfei • x) 
-3A 2 (2|A| 2 + \A 2 \ 2 + 2\A 3 \ 2 ) exp(ik 2 ■ x) 
-3A 3 (2|A| 2 + 2\A 2 \ 2 + \A 3 \ 2 ) exp{ik 3 • x) 
-QA 2 A* 3 ip exp(ikx ■ x) — 6AlAlip exp(ik 2 ■ x) 

— QA\A* 2 ip exp(ik 3 ■ x). 



(13) 



Aligning the secular terms in Eq. along each basis, and scaling back to the original 
variables, we obtain the amplitude equations as 



Of 



CiA - 3A 1 [lAi] 2 + 2\A 2 \ 2 + 2\A 3 \ 2 ) - QipA* 2 A* 3 . 



together with the appropriate permutations for A 2 and A3. 
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(a) 
t = 56.04 



(b)t = 
184.04 



(c)t = 
720.04 



(d) 
t = 56.04 



(e)t = 
184.04 



(f)t = 
720.04 



FIG. 1: Comparison of heterogeneous nucleation and growth in the PFC equation Eq. © (panels 
(a)-(c)) and its RG-generated mesoscale counterpart Eq. @3J) (panels (d)-(f)). The order parameter 
is shown at the times indicated starting from the same initial condition with ip = 0.285 and 
r = -0.25. 



Fig. (0) shows the time evolution for the nucleation and growth of a two-dimensional film, 
calculated using the PFC equation and its RG-generated mesoscale counterpart. Starting 
from the same initial condition of randomly-oriented seeds, crystalline domains grow, collid- 
ing to form a polycrystalline microstructure. The solutions from the two different computa- 
tional algorithms are essentially indistinguishable. Without the PFC formulation, it would 
not have been possible to capture successfully the formation of a polycrystalline material, 
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with grain boundaries and other defects. We conclude that it is indeed possible to compute 
large scale microstructure from effective equations at the mesoscale. 

IV. COMPUTATIONAL EFFICIENCY 

In order for our scheme to be computationally efficient, we need to establish that the 
solutions for the fundamental mesoscale variables are indeed slowly varying. In Fig. (J2J) are 
shown the amplitude and phase gradient of one of the components during the computation 
of the two-dimensional grain growth. It is evident that the variables are indeed essentially 
uniform, except near the edges of the grains. 




(a) (b) 

FIG. 2: Color map of the spatial variation of the amplitude (a) and phase gradient (b) of the 
solution displayed in Fig. (JJ) at time t = 720.04. 

In order to exploit this property computationally, we need to work with a formulation 
which is independent of the particular orientation of our reference directions. To this end, 
we reformulate the RG-PFC equations for the complex amplitudes Aj etc. in terms of their 
real amplitude and phase variables, denoted by tyj > and $j . Expanding out the terms in 
Eq. ()14j) . and equating real and imaginary parts, we obtain equations of motion for tyj and 
V$j which can readily be solved by adaptive mesh refinement, to be reported elsewhere. 
This formulation is also important in treating the beats that can arise if the crystallographic 
axes of a grain are not collinear to the basis axes used in the numerical solution. Such beats 
can be dealt with by either adaptive mesh refinement or the polar coordinate formulation, 
and will be discussed in detail elsewhere. As we have previously shown |4], adaptive mesh 
refinement algorithms scale optimally, with the number of computations being proportional 



S 



to the number of mesh points at the finest scale, i.e. to the grain boundary length (two 
dimensions) or surface area (three dimensions). 

In summary, we have shown that multiscale modeling of complex polycrystalline materials 
microstructure is possible using a combination of continuum modeling at the nanoscale using 
the PFC model, RG and related techniques from spatially-extended dynamical systems 
theory. The PFC model can be extended to include three dimensions, multi-component 
systems, thermal fields and realistic atomic correlations. Our analysis is ideally-suited for 
efficient adaptive mesh refinement, thus enabling realistic modeling of large-scale materials 
processing and behavior. 
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